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Abstract 

The sum of correlated gamma random variables appears in the analysis of many wireless commu- 
nications systems, e.g. in systems under Nakagami-m fading. In this Letter we obtain exact expressions 
for the probability density function (PDF) and the cumulative distribution function (CDF) of the sum 
of arbitrarily correlated gamma variables in terms of certain Lauricella functions. 

Index Terms 

Gamma variates, Nakagami-m fading. Outage Probability, Lauricella Functions. 

I. Introduction 

Many of the performance analysis problems in the scope of wireless communications theory 
require determination of the statistics of the sum of the squared envelopes of Nakagami-m faded 
signals or, equivalently, the sum of gamma random variables since the square of a Nakagami-m 
random variable follows a gamma distribution [1]. 

Some expressions are available in literature for the probability density function (PDF) of the 
sum of gamma random variables, e.g. see [1] and the references cited herein. These expressions 
are frequently in the form of infinite series, including the general expression for arbitrary 
correlation provided in [1]. In this Letter we revisit the result derived in [1]; we show that, 
under the same assumptions, the PDF and the cumulative distribution function (CDF) of the 
sum of correlated gamma random variables can be expressed in a compact form by certain 
Lauricella functions. To the best of the author's knowledge, the expressions obtained here are 
novel. 

*This work is partially supported by the Spanish Government under project TEC2007-67289/TCM and by AT4 wireless. 
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II. Analytical Results 

For clarity we will use the notation adopted in [1]. We say X follows a gamma distribution 
with parameters a > and /3 > if the PDF of X is given by 

a-l -x/f3 

where r( ) is the gamma function and C/( ) is the unit step function. The shorthand notation 
X ~ Q{a,^) will be used to denote that X is gamma distributed with parameters a and ^. 

The key idea of this Letter is the following. In [1] the authors extended the Moschopoulos' 
Theorem [2]. Interestingly, they observed the similarity between the MGF of the correlated case 
and the independent case. Then, the Moschopoulos technique of inverting the MGF was again 
adopted for the correlated case. Here, we start from the MGF of the correlated case, but instead 
of using the Moschopoulos technique, we establish a connection between the MGF and certain 
Lauricella functions. The first who established a connection of this type was Kabe [3], within 
the context of independent gamma random variables. The mathematically precise statements are 
given below. 

Lemma 1: Let {Xn}n=i be a set of N correlated not necessary identically distributed gamma 
random variables with parameters a and respectively, [i.e., X^ ~ G{a, Pn)] and let pij denote 
the correlation coefficient between Xj, and Xj, i.e., 

cov (Xi, Xj) 

Vvar(^i)var(Xj) (2) 

then the CDF of r = YZ=i 

can be expressed as 



Na 

Fy {y) ^ 



det(^)"r (1 + Na) 



(3) 



where $2^^ is the confluent Lauricella function [4][5], and {Xn}n=i are the eigenvalues of the 
matrix A — DC where D is the N x N diagonal matrix with the entries {/3n}n=i and C is the 



January 13, 2013 



DRAFT 



3 



N X N positive definite matrix defined by 

/ 1 



1 



\ VPni 



The PDF of Y is given by 



friy) 



-1+Na 



dct(A)°r (Na) 



^/PlN ^ 
y/P2N 

1 J 



(4) 



(5) 



Proof: See Appendix I. 



The expressions derived in Lemma 1 are compact and can be frequently reduced to simpler forms 
using the properties of the function $2^\ In particular, since $2^^ = iFi, [i.e., equivalent to the 
confluent hypergeometric function] one can check that for a = m, /3i = 7/m and A?^ = 1 the 
expressions derived here reduce to the well-known CDF and PDF of the square of a Nakagami-m 
random variable. For reduction formulas, integral representations and integrals involving $2^'*, 
the reader should refer to [4]. Note that the CDF expression given Lemma 1 allows us to compute 
the outage probability of maximal ratio combining (MRC) over correlated Nakagami-m fading 
channels. 



in. Conclusions 

In this Letter, compact expressions have been derived for the sum of arbitrarily correlated 
gamma random variables. Such expressions have both theoretical and practical value, and are 
applicable in a vast range of wireless communications problems. 

Appendix A 
Proof of Lemma 1 

For an arbitrary function (f){x) we denote the Laplace transform as C[(}){x); s]. As in [1], we 
define the MGF of Y as My{s) = E[e'y] = jC[fY{y)] s]. We know that the MGF of Y is given 
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N 



MY{s) = l[{l-XnS)-''. (6) 



n=l 



Therefore, we can write the CDF of Y in the following form 

£[FY{y);s] = -C[fY{y);s] 

S 




(7) 



s 

Then, after identifying (7) with [5, p. 222, eq. 5], the CDF of Y is obtained. To derive the 

expression for the PDF we can write 



N 



n=\ 

N , ^Oi\ ( f N \ 



(8) 



s 



n=l 



-a 



\n=l 
n=l \ 

Again, after identifying this expression with [5, p. 222, eq. 5], the PDF of Y is obtained. 



References 

[1] M.-S. Alouini, A. Abdi and, M. Kaveh, "Sum of gamma variates and performance of wireless sommunications systems 
over Nakagami fading channels ," IEEE Trans. Veh. TechnoL, vol. 50, pp. 1471-1480, Nov. 2001. 

[2] P. G. Moschopoulos, "The distribution of the simi of independent gamma random variables ," Ann. Inst. Statist. Math. (Part 
A), vol. 37, pp. 541-544, 1985. 

[3] D. G. Kabe, "On the exact distribution of a class of multivariate test criteria," Ann. Math. Statist., vol. 33, pp. 1197-1200, 
1962. 



January 13, 2013 



DRAFT 



5 

[4] H. Exton, Multiple Hypergeometric Functions and Applications, (Halsted Press) John Wiley & Sons, 1976. 
[5] A. Erdelyi, Tables of Integrals Transforms, vol. I, McGraw Hill, New York, 1954. 



January 13, 2013 



DRAFT 



